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Preface 


The  idea  behind  this  thesis  came  from  articles  printed 
in  the  bulletin  of  the  Space  Studies  Institute  of  Princeton, 

New  Jersey.  The  members  of  SSI  have  been  most  helpful  in 
directing  me  to  sources  of  information  on  the  External  Tank 
and  its  potential  use  in  space. 

In  order  to  perform  an  orbital  analysis  a  method  of 
determining  a  satellite’s  position  at  some  future  time  is 
required.  A  sophisticated  orbit  prediction  computer  program 
called  the  Artificial  Satellite  Analysis  Program  (ASAP) 
fulfills  this  need.  I  am  grateful  to  my  advisor,  Captain 
Rodney  Bain,  for  instructing  me  in  its  operation.  His 
enthusiasm  in  the  subjects  of  astrodynamics  and  celestial 
mechanics  provided  a  key  motivation  for  this  work. 

As  section  leader  for  GA-87D/88M,  I  would  like  to  thank 
the  other  members  for  their  help  during  our  tour  at  AFIT  and 
for  their  assistance  in  getting  me  up  to  speed  after  my 
twelve  year  academic  layoff.  And  finally,  I  wish  to  express 
my  love  and  respect  to  my  wife,  Gail,  who  has  given  steadfast 
support  and  encouragement  during  these  challenging  scholastic 
pursuits. 
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Abstract 

By  flying  a  different  launch  profile  ,  it  is  possible 
for  the  Space  Transportation  System's  Orbiter  to  bring  the 
External  Tank  directly  into  space.  Many  studies  by  NASA  and 
priv'ate  industry  have  detailed  the  potential  on-orbit  uses  of 
an  External  Tank.  However,  at  Space  Shuttle  operating 
altitudes,  an  orbiting  tank  will  experience  multiple  environ¬ 
mental  forces  resulting  in  its  decay  into  the  lower  atmos¬ 
phere  and  eventual  re-entry. 

This  thesis  conducts  a  preliminary  study  of  a  single 
External  Tank  in  low  Earth  orbit.  Criteria  for  a  parking 
orbit  are  defined  and,  using  an  orbit  prediction  computer 
program  with  atmospheric  drag  and  gravitational  perturbations 
included,  a  search  is  made  for  the  lowest  initial  altitude 
that  will  allow  the  External  Tank  to  remain  in  this  orbit 
window.  The  starting  altitude  that  meets  the  orbit  require¬ 
ments  is  found  to  be  within  reach  of  the  Shuttle’s  capabili¬ 
ties.  The  orbital  elements  of  this  parking  orbit  are  then 
analyzed  and  a  method  for  quick  calculation  of  these 
parameters  is  devised.  An  evaluation  of  the  factors  that 
affect  the  orbital  contraction  of  an  External  Tank  is  also 
performed.  The  atmospheric  density  and  the  tank  character¬ 
istics  can  both  contribute  to  high  orbital  decay  rates. 


ORBITAL  ANALYSIS  OF  A  STS  EXTERNAL  TANK 


IN  LOW  EARTH  ORBIT 


I .  INTRODUCTION 


One  of  the  main  goals  of  the  Space  Transportation  System 
(STS)  is  to  reduce  the  cost  of  delivering  payloads  into 
orbit.  Reuse  of  the  Solid  Rocket  Boosters  and  the  Orbiter 
itself  is  a  major  factor  in  achieving  this  goal  but  the 
original  Space  Shuttle  designers  elected  to  make  the  External 
Tank  ( ET )  an  expendable  element  of  the  STS.  Hence,  on  each 
flight,  the  ET  is  jettisoned  by  the  Orbiter  for  a  controlled 
entry  into  the  Indian  or  Pacific  Ocean.  However,  this 
component  of  the  STS  does  not  have  to  go  to  waste.  When  cast 
off,  the  External  Tank  has  98%  of  the  energy  needed  to  insert 
it  into  orbit.  A  more  efficient  launch  trajectory,  called 
direct  injection,  could  be  flown  where  the  Space  Shuttle  Main 
Engines  would  boost  the  Orbiter  with  its  attached  External 
Tank  into  a  standard  Shuttle  orbit  (See  Figure  1).  The 
orbiter  would  then  leave  the  ET  and  continue  on  its  scheduled 
mission  ( 8 : 1 , I  -  1  ) . 

Why  bring  the  External  Tank  into  orbit?  Many  on-orbit 


applications  of  an  ET  have  been  proposed  by  government  and 
private  industry  groups.  This  paragraph  presents  a  review  of 


Figure  1.  Direct  Injection  Trajectory 

some  of  these  proposals  as  discussed  in  Reference  8.  The 
External  Tank  has  approximately  53,000  lbs  of  aerospace  grade 
aluminum.  The  tank  could  be  disassembled  in  orbit  with  the 
pieces  used  to  construct  large  space  structures.  It  is  also 
possible  to  melt  the  tank  aluminum  for  on-orbit  manufacturing 
uses,  or  leave  the  ET  whole  and  use  it  as  a  basis  for  a 
space  station.  Each  ET  consists  of  a  liquid  hydrogen  tank 
plus  a  smaller  liquid  oxygen  tank;  both  capable  of  on-orbit 
storage  of  cyrogenic  fuels  and  other  volatiles.  Studies  have 
been  done  to  investigate  the  employment  of  tethers  in  space. 
The  potential  uses  of  the  ET  connected  to  a  tether  range  from 


2 


momentum  exchange  with  the  Shuttle  to  electical  generation 
with  a  conducting  tether.  However,  prior  to  utilizing  an 
External  Tank  in  space,  a  major  obstacle  must  be  addressed; 
how  to  prevent  it  from  decaying  out  of  orbit  prematurely. 

All  objects  in  low  Earth  orbit  (LEO)  experience 
atmospheric  drag.  Work  is  required  to  push  the  air  molecules 
out  of  the  way  which  reduces  the  kinetic  energy  of  a 
satellite.  This  effort  causes  the  orbit  to  shrink  where  the 
orbital  velocity  requirement  is  greater.  This  increase  in 
speed  combined  with  the  higher  density  of  the  lower 
atmosphere  results  in  an  increased  drag  force.  And  the  cycle 
continues  until  the  satellite  falls  out  of  orbit  (25:296). 
This  phenomona  was  observed  with  the  Skylab  space  station. 
Skylab’s  orbit  decayed  due  to  increased  solar  activity  that 
affected  the  density  of  the  atmosphere.  The  additional  drag 
on  the  spacecraft  led  to  its  earlier  than  planned  re-entry 
(24:39).  An  orbiting  External  Tank  would  suffer  the  same 
fate.  An  unexpected  re-entry  would  defeat  the  purpose  for 
bringing  the  ET  into  space  and  produce  the  possibility  of 
raining  large  pieces  of  a  disintegrating  tank  over  populated 
areas.  A  method  of  orbit  maintenance  is  desirable  but  the 
first  tanks  taken  into  orbit  may  not  be  equipped  with  small 
thrusters  to  accomplish  this  task.  And  possibly,  ET 
end-users  would  make  use  of  the  tank  before  it  began  a  final, 
blazing  plunge  through  the  lower  layers  of  the  atmosphere. 

The  thrust  of  this  thesis  is  to  determine  a  parking 
orbit  for  a  single  External  Tank,  awaiting  its  future 
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purpose.  Two  qualifying  assumptions  are  made  to  form  a  basis 
for  this  investigation.  First,  the  time  horizon  chosen  is  90 
days.  This  three  month  window  is  selected  to  represent  the 
amount  of  time  allowed  for  the  ET  to  remain  in  its  parking 
orbit.  It  is  also  assumed  that  the  maximum  tolerable 
altitude  loss  by  the  ET  over  the  90-day  period  would  be  25 
kilometers.  These  two  assumptions  thus  produce  a  platform 
from  which  to  answer  the  primary  question:  What  is  the 
lowest  initial  altitude  for  deploying  the  External  Tank 
sufficient  to  meet  this  parking  orbit  requirement? 

An  orbit  prediction  computer  program,  coupling  a 
gravitational  field  and  a  simple  atmosphere  model,  is 
employed  for  this  External  Tank  orbital  analysis.  Typical 
Shuttle  orbit  parameters  of  28  inclination  and  zero 
eccentricity  (circular  orbit)  are  used  in  the  search  for  a 
minimum  parking  orbit.  The  Keplerian  elements  a,e,i,Q  and 
the  geocentric  altitude  of  this  orbit  are  then  analyzed  over 
the  90-day  period  and  a  simple  method  is  devised  allowing 
quick  calculation  of  these  ET  orbital  parameters.  Finally, 
various  aspects  of  the  orbit,  Earth’s  gravitational  field  and 
atmosphere,  and  External  Tank  properties  are  investigated  to 
determine  their  relationship  with  orbital  decay.  Chapter  V 
details  the  results  of  this  examination  of  an  External  Tank 
in  low  Earth  orbit. 


I I .  COMPUTER  PROGRAMS 


Several  computer  pro>?rams  were  used  in  formulating  and 
analyzing  the  External  Tank  orbital  data.  The  data  itself 
was  generated  by  the  Artificial  Satellite  Analysis  Program 
(ASAP)  from  Reference  16  while  data  analysis  and  graphing  was 
accomplished  using  various  programs  and  software  packages. 


Artificial  Satel 1 i te  Anal ys i s  Program  ( ASAP ) 

ASAP  is  a  general  orbit  prediction  program  written  in 
the  standard  ANSI  77  FORTRAN  language.  Though  originally 
configured  to  run  on  an  IBM  microcomputer,  ASAP  was  converted 
to  operate  on  AFIT’s  mainframe  systems.  This  program  uses 
Cowell’s  method  of  special  perturbations  in  formulating  the 
differential  equations  of  motion  ( EOM ) . 

Cowell’s  method  works  with  cartesian  coordinates  (x,y,z) 
and  their  time  derivatives  (x,y,z).  The  EOM  can  be  expressed 


y  =  V  y  = 


z  r  Vi 


p  -  +  perturbations 

r 


-  p  •  ^  +  perturbations 
r 


-  p  -  +  perturbations 


where 


r 


(x^  + 


2,1/2 

Z  ) 


M  =  G  M 

ASAP’s  numerical  integrator  requires  the  equations  of 
motion  be  written  as  a  system  of  6  first  order  differential 
equations.  The  integrator  is  an  8th  order  Runge-Kutta  single 
step  integration  scheme  with  variable  step-size  control.  The 
perturbations  handled  by  ASAP  include  the  primary  body 
disturbing  function,  atmospheric  drag  and  solar/lunar 
effects.  For  the  sake  of  completeness,  a  general  discussion 
of  the  gravitational  and  atmospheric  drag  perturbations  is 
provided  in  Chapter  III.  The  method  of  their  implementation 
in  ASAP  is  then  covered  in  Chapter  IV. 

In  addition  to  basic  constants  for  the  central  planet 
the  following  are  user  inputs  to  ASAP: 


1)  Coefficients  of  Spherical  Harmonics 

2)  Atmospheric  Data 

A)  Reference  Altitude  (km) 

B)  Scale  Height  at  Reference  Altitude  (^n>) 

C)  Density  at  Reference  Altitude  (kg/km  ) 

3)  Satellite  Data 

2 

A)  Area, normal  to  flow  (km  ) 

B)  Mass  (kg) 

C)  Coefficient  of  Drag 

4)  Initial  Orbital  Elements 

A)  Semi-Major  Axis  (km) 

B)  Eccentricity 

C)  Inclination  (deg) 

D)  Longitude  of  Ascending  Node  (deg) 

E)  Argument  of  Periapsis  (deg) 

F)  Mean  Anomaly  (deg) 
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5)  Starting  Day  and  Time 


6)  Stopping  Day  and  Time 

7 )  Time  Step 


At  each  time  step  ASAP  integrates  the  equations  of 
motion,  Eqs  (1),  and  calculates  the  Keplerian  elements. 

Since  perturbations  have  been  included  in  the  EOM ,  the 
orbital  elements  will  not  remain  constant  but  be  functions  of 
time.  The  orbit  described  by  these  elements  is  the  two  body 
orbit  a  satellite  would  follow  if  all  the  perturbations  were 
removed  at  that  moment  in  time  (28:29-30).  These  instanta¬ 
neous  values  of  the  orbital  elements  are  used  as  a  basis  for 
calculating  the  satellite’s  orbit  parameters  at  the  next  time 
period.  At  each  step  ASAP  generates  the  following  output: 


1  )  Day 

2 )  Hour 

3)  Semi-Major  Axis  (km) 

4 )  Eccentricity 

5)  Inclination  (deg) 

6)  Longitude  of  Ascending  Node  (deg) 

7)  Argument  of  Periapsis  (deg) 

8)  Mean  Anomaly  (deg) 

9)  Geocentric  Altitude  (km) 


Other  Prog rams 

The  Artificial  Satellite  Analysis  Program,  as  it 
calculates  the  orbital  elements,  produces  an  ASCII  file 
which  can  be  imported  to  spreadsheet  type  software.  The 
database  management  capabilities  of  this  program  allows  rapid 
analysis  of  the  orbital  parameters.  For  the  External  Tank 
orbit  study  a  two  hour  time  step  was  employed  which  produced 


12  orbit  element  calculations  per  day  for  the  three  month 
period.  The  value  of  the  parameters  at  the  zero  hour  of  each 
day  was  extracted  for  integration  into  the  graphics  portion 
of  the  software.  The  orbital  element  time  histories  were 
produced  with  this  data.  A  linear  regression  function  was 
another  program  feature  used.  As  will  be  discussed  in 
Chapter  V,  this  capability  provided  a  conv'enient  method  for 
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III.  ORBIT  PERTURBATIONS 


A  satellite,  orbiting  a  spherical  planet  of  homogeneous 
structure  and  with  no  outside  environmental  disturbances, 
would  continue  forever  in  an  idealized  Keplerian  orbit;  its 
orbital  elements  remaining  constant  at  their  initial  values. 
In  real  celestial  systems  however,  there  are  interfering 
influences  capable  of  perturbing  the  satellite’s  motion: 


1  ) 

primary 
2  ) 

3  ) 

4  ) 

5  ) 


the  irregular  shape  and  variable  composition  of  the 

the  atmosphere  of  the  primary 
the  magnetic  field  of  the  primary 
the  solar  radiation  pressure 

the  gravitational  fields  of  other  nearby  bodies 


This  thesis  will  concern  itself  primarily  with  the  first  two 
perturbations  listed. 

The  magnitude  of  a  perturbing  force  gives  no  hint  as  to 
its  effect  on  a  specific  orbital  element  (9:973).  For 
instance,  the  Earth’s  gravitational  attraction  is  on  the 
order  of  one  dyne,  but  because  it  is  a  conservative  force, 
produces  no  appreciable  long  term  perturbation  upon  a 
satellite’s  semi-major  axis.  Conversely,  the  atmospheric 
drag  (on  the  order  of  10  ^  dyne)  acting  on  a  low  orbit 
satellite  causes  a  noticeable  secular  decrease  in  the 
semi-major  axis. 

The  next  two  sections  will  discuss  in  detail  the  two 
principal  perturbations  experienced  by  a  satellite  in  low 
Earth  orbi t . 
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Gravi  tation a  1  r bat  ions 

As  experienced  by  a  satellite  in  low  orbit,  Earth’s 
gravitational  field  does  not  exhibit  spherical  symmetry.  The 
equatorial  bulge,  continent  land  masses  and  irregular  ocean 
floors  all  contribute  to  Earth’s  deviation  from  a  Newtonian 
point  mass  and  produce  a  non-spherical  shape  to  Earth’s 
gravitational  field  (28:48).  As  a  satellite  circles  the 
globe  it  encounters  a  variable  force  of  gravity  leading  to 
changes  in  its  orbital  parameters.  A  formulation  is  needed 
to  model  this  perturbing  effect. 

The  attractive  force  between  two  bodies  of  mass  M  and  m, 
separated  by  distance  r,  is  governed  by  Newton’s  universal 
law  of  gravity; 

M  m 

F  =  G  -  (  2  ) 

2 

r 

Combining  Eq  (2)  with  Newton’s  second  law 

F  =  m  a  (  3  ) 

provides  the  acceleration  of  body  ro  with  respect  to  the 
two-body  system’s  center  of  mass: 

G  M 

a  =  -  (  4  ) 

2 

r 

For  the  case  where  m  <<  M,  the  two-body  system’s  center  of 
mass  can  be  conveniently  considered  at  the  center  of  body  M. 
Instead  of  working  with  an  acceleration  in  determining  the 
gravitational  disturbance,  it  is  simpler  to  convert  the 


acceleration  to  a  vector  and  express  this  vector  as  a  poten¬ 
tial  function  (14:1-2). 


The  vector  a  is  obtained  by  represent  ng  the  accel¬ 
eration  as  a  gradient  of  a  scalar,  V,  defined  as  the 
potential  (per  unit  mass): 

a  =  ^  V  (  5  ) 


where 


r 

This  acceleration  a  is  used  in  developing 
fundamental  partial  differential  equation 
fields  (6:277,279): 


(  6  ) 

Po i sson  ’  s 

for  gravitational 


V^V  =  -4nGp 


(  7  ) 


where  p  =  density  of  M.  For  a  satellite  orbit  outside  the 
body  M,  the  density  equals  zero  which  leads  to  Laplace's 
equal  ion: 

17^  V  =  0  (8) 


In  rectangular  coordinates,  Laplace’s  equation  becomes 


17^  V  =  G  M 


2  ?  2 
+  3(x  1  y  t  z  ) 

3  5 

r  r 


0 


(  9  ) 


Since  the  body  M  (Earth,  in  this  case)  is  basically  spherical 

in  shape,  transforming  Eq  (9)  into  spherical  coordinates 

will  aid  in  the  derivation.  This  can  be  accomplished  with 

X  =  r  cos<P  cosX 

y  =  r  cos't'  cosX 

z  -  r  s in<t>  (10) 


where 


r  ^  radial  distance  from  coordinate  frame  origin 
<t>  ^  latitude 
X  =  longitude 

Thus,  Laplace’s  equation  transforms  to  (13:2-3) 


o  2 

r  = 


dr 


2  dV 


dr 


1 


cos4>  d4> 


cos4> 


d\^ 


d<t> 


1 

cos^<P  dX^ 


(11) 


Any  solution  V  to  Eq  (11)  is  referred  to  in  literature  as  a 
spherical  harmonic  but  a  closed  form  solution  does  not  exist. 
Instead,  a  general  infinite  series  expansion  for  V  can  be 
derived  (14:4-6): 


f  1 


- 1  - 1 


V  =  y  y  r  ‘  ‘  P  (sin<^’)  IC  cus  mX  +  S  sin  mX 

1  ni  I  1  m  1  tn 

1=0  m=0 


Eq  (12)  involves  an  associated  Legendre  function 

k 

P  (sin<l>)  =  cos”'<P  y  T  sin* 

1  m  1  m  I 

t=0 

where  k  =  integer  part  of  (l-m)/2  and 


(12) 


(  13  ) 


( -1  )  ( 2  1-2  t  )  : 


2 '  t !  (  i-t )  :  (  1  -m-2 1 )  : 


(  14  ) 


The  constants  C  and  S  characterize  the  mass  distribution 

1  m  1  fn 

of  body  M  and  the  terms  1  and  m  are  the  order  and  degree  of 
the  potential.  In  this  thesis  the  notation  3x3  refers  to  a 
potential  of  order  3  and  degree  3.  To  represent  a  model  of 
the  Earth’s  geopotential,  Eq  (12)  can  be  expressed  as  (6:284) 


*•*  1 

V  =  -~  y  y  P  (sin<l>)  fc  cos  mX  +  S  sin  mX]  (15) 

r  ^  ^  -  1  m  1  m  Im  ^ 


I «0  msO  r 


By  including  the  Earth's  gravitational  parameter  ^J  and 
equatorial  radius  R  ,  C  and  S  are  redefined  from  Eq  (12) 

e  1  m  1  m 

to  make  them  dimensionless. 

When  m  =  0  Eq  (15)  simplifies  to 


V  =  ^  I 


P,  (sin<i>)  (C,J 


The  geopotential  components  of  this  equation  are  referred  to 
as  zonal  harmonics  and  are  due  to  variations  in  meridian 
ellipticity.  When  m  =  1,  the  components  of  Eq  (15)  are 
called  sectoral  harmonics  and  are  caused  by  longitudinal 
variations  in  the  shape  of  the  Earth.  Tesseral  harmonics  are 
concerned  with  cases  where  m  <  1  (6:271).  (Note:  For 

all  mil,  the  associated  Legendre  function  P  (sin<J>)  equates 

1  rt) 

to  zero.)  Values  of  the  potential  coefficients  C  and  S 

i  m  i  m 

for  Earth  were  determined  from  repeated  observations  of 
various  satellite  orbits.  By  studying  the  long  term 
perturbations,  the  coefficients  of  the  zonal  harmonics  have 
been  found  and  the  sectoral  and  tesseral  l..-.’'monic 
coefficients  were  discovered  from  analysis  of  short  term 
orbit  perturbations.  These  geopotential  constants  have  been 
calculated  through  many  orders  and  degrees. 

The  3x3  geopotential  coefficients  of  concern  to  this 
thesis  are  listed  in  Table  I  (6:285). 


TABLE  1 


Coefficients  of  Geopotential  Harmonics 
>6 


m 


C.  X  10 

1  nt 


S  X  10 

1  in 


2  -1082.645  2  2  2.379  -1.351 

3  2.546  3  1  1.936  0.266 

3  2  0.734  -0.538 

3  3  0.561  1.620 


When  both  1  and  m  equal  zero,  Eq  (15)  reduces  to  the  basic 
potential  for  a  spherically  symmetric  Earth: 


V  =  -y_  =  -G  M 

r  r 

meaning  1  and  0  (28:56).  By  placing  the 

coordinate  frame  at  the  center  of  th  Earth,  coefficients  C 

C,  .  will  be  zero  and,  from  empirical  data,  C  and 
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have  been  shown  to  be  exceedingly  small  (6:285).  The 
sectoral  and  tesseral  harmonic  coefficients  C  and  S  in 

1  m  1  m 

Table  I  are  related  to  the  contants  C  and  S  of  Eq  (15) 

1  m  1  m 


(6) 


lo’ 


by 
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I  m 


T  C 


1  m 


S  = 

1  m 


T  S 


I  m 


(17) 


where 


T  = 


(214-1  )  (  l-m) 
(  1+m)  ! 


1/2 


(  18  ) 


At  this  stage  the  geopotential  has  been  developed  into  a 
form  used  by  the  Artificial  Satellite  Analysis  Program.  As 
it  is  of  interest  to  this  thesis  to  examine  the  long  period 


perturbating  effects  of  Earth’s  gravitational  field,  the  time 
histories  of  the  External  Tank  orbital  elements  are  generated 


by  ASAP  and  then  a  curve  fitting  to  this  data  is  performed. 
This  approach  to  determine  the  long  term  trends  is  a  reverse 
procedure  from  an  analytical  technique  called  the  method  of 
averages  (12:1).  In  this  process  the  disturbing  function, 
i.e.,  the  geopotential,  is  averaged  prior  to  its  use  in  the 
equations  of  motion.  The  is  accomplished  by  using 


where  the  potential  is  averaged  over  one  orbit  period  (13:3). 
However,  as  Eq  (15)  now  stands,  the  geopotential  is  a 
function  of  radius  r,  latitude  <t>  and  longitude  X.  For  use  in 
the  method  of  averages  it  is  necessary  to  convert  Eq  (15)  to 
the  Keplerian  elements  a,e,i,0,(J  and  M.  This  conversion  is 
presented  in  Appendix  A. 


Atmospheric  Drag  Perturbations 

As  a  satellite  passes  through  the  rarified  regions  of 
the  upper  atmosphere  it  is  subjected  to  an  aerodynamic  force 
which  perturbs  its  orbital  motion.  This  force  can  be 
separated  into  two  quantities:  (1)  the  component  opposite 
the  direction  of  motion  called  drag,  (2)  the  component 
perpendicular  to  the  satellite’s  flight  path.  This  second 
quantity  generally  docs  not  pass  through  the  satellite’s 


center  of  mass  and  can  subsequently  be  divided  into  a  lifting 
force  and  a  turning  moment  about  the  center  of  mass.  Accord¬ 
ing  to  Reference  15  the  aerodynamic  lift  can  be  neglected  in 
an  orbital  analysis  as  the  lift  to  drag  ratio  for  satellites 
is  very  small  (<0.1).  Also, if  uncontrollable,  a  satellite 
would  rotate  the  lift  vector  through  all  possible  directions 
over  a  period  of  time  thus  producing  a  resultant  force  of 
zero.  These  assumptions  for  ignoring  the  lift  force  are 
particularly  reasonable  for  cylindrical  satellites  with 
length  to  diameter  ratios  greater  than  1  (15:13).  (The  ratio 

of  1 /d  for  the  External  Tank  is  5.6.)  In  addition,  a 
cylindrical  satellite  in  a  gravity  gradient  orientation  with 
its  longitudinal  axis  pointing  at  the  Earth’s  center  would 
not  generate  any  lift. 

In  aerodynamics  the  conventional  drag  equation  is 

D  =  1  /2  p  V^Cd  A  (  20  ) 

where 

D  ^  drag  force 
p  ?  atmospheric  density 

V  ^  satellite  velocity  relative  to  the  atmosphere 
('d  drag  coefficient 

A  ^  satellite  area  normal  to  air  flow 


This  formulation  of  the  drag  equation  can  alternately  be 
expressed  as  a  drag  deceleration: 


The  ba 1 1 i st i c 
combining  the 


D 


Cd  A 


m 


2m 


p  V 


B  p  V 


(  21  ) 


coefficient  R  is  a  convenient  method  of 
satellite  parameters  ( A , m , Cd  )  into  one  term  for 


parameteric'  studies  (1  1:183-184). 
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The  elements  of  the  drag  equation  will  now  be  examined 
with  respect  to  an  Earth  orbiting  spacecraft. 


Atmospheric  Dens i ty .  The  major  property  of  the  Earth’s 
atmosphere  important  to  low  orbit  satellites  is  its  density. 
The  density  in  the  upper  regions  of  the  atmosphere  is  an 
exceedingly  variable  quantity  with  solar  activity  being  the 
main  cause  of  this  irregularity.  Observations  of  the  orbits 
of  early  satellites  were  used  in  analyzing  these  atmospheric 
density  fluctuations.  (9:978-979). 

One  basic  air  density  difference  exists  between  the 
sunlit  and  nighttime  sides  of  the  Earth.  See  Figure  2  from 
Reference  11.  This  hump  of  density,  the  diurnal  bulge, 


Figure  2.  Diurnal  Atmospheric  Bulge 


occurs  about  two  hours  after  midday.  The  change  in  density 
between  day  and  night  is  small  for  altitudes  of  250  km  and 
below,  however,  above  this  altitude  large  variations  can  take 
place.  For  example,  at  600  km  ,  the  maximum  daytime  density 
is  about  8  times  greater  than  the  nighttime  minimum  (15:18). 

Along  with  this  daily  density  variation  there  are  four 
other  solar  influences.  First,  there  exists  an  irregular 
day-to-day  deviation  in  density  due  to  ephemeral  solar 
disturbances  such  as  short  term  flares.  A  second  density 
variation  occurs  with  a  27-day  cycle.  This  effect  is  due  to 
the  axial  rotation  of  the  Sun  with  respect  to  the  Earth, 

Next,  the  properties  of  the  atmosphere  also  respond  to  the 
10-11  year  sunspot  cycle.  And  fourth,  a  seasonal  oscillation 
in  density  appears  with  a  minimum  occurring  during  July  and 
January  and  a  maximum  appearing  in  October  and  April.  These 
seasonal  effects  are  largely  due  to  the  Earth’s  orbit  around 
the  Sun  ( 15 : 18-19  )  . 

Through  many  orbital  studies  of  Earth  satellites  the 
atmosphere  and  its  properties  have  been  analyzed  in  much 
detail.  From  these  observations  the  atmospheric  density  has 
shown  to  vary  exponentially  (as  a  first  approximation)  with 
height  (15:20,22).  Figure  3  is  a  logarithmic  plot  of  air 
density  versus  altitude  that  demonstrates  this  effect 
(21:20).  The  relatively  straight  line  relationship  between 
200  and  800  km  allows  the  atmospheric  density  to  be  expressed 
by  a  simple  exponential  law  (15:22): 

p  =  exp  I  -(h-h„)/H  )  (22) 


Figure  3.  Atmospheric  Density  Variation  with  Altitude 


where 

h  =  altitude 

h  =  reference  altitude 
0 

p  H  density  at  reference  altitude 
0 

H  =  scale  height 

The  scale  height  is  the  vertical  distance  in  which  the  air 
density  changes  by  a  factor  of  exponential  e  (22:233).  The 
value  of  H  is  dependent  on  altitude,  increasing  slowly  with 
height.  H  can  be  derived  in  the  following  manner  (27:4). 
Two  basic  equations  play  a  part  in  determining  air 


density  as  a  function  of  altitude. 


First,  the  perfect  gas 


VS-Vi 


law  relates  the  atmospheric  properties  pressure  P,  density  p 


and  temperature  T: 


where 


P  =  p 


R  =  universal  gas  constant 

M  =  mean  molecular  weight  of  the  atmosphere 
The  second  expression  is  the  hydrostatic  equation  which 
relates  the  rate  of  change  in  pressure  to  the  increased 
weight  of  the  supported  atmosphere: 

dP  =  -  p  g  dr  (2 

where  g  s  gravitational  acceleration.  Logarithmic  differ¬ 
entials  of  Eq  (23)  give 


When  combined  with  Eq  (24)  this  produces 


-  ^  dr 


where  |3  is  the  terms  in  the  brackets  of  Eq  (26)  and 
represents  the  inverse  of  the  scale  height.  Integrating 
Eq  (27)  and  substituting  H  for  1/P  yields  the  density 
exponential  equation  Eq  (22). 

Sever a  1  different  density  profiles  can  be  considered 
depending  on  assumptions  on  P  (27:4-5): 

1)  Strictly  exponential  atmosphere:  P  is  assumed 
constant  throughout  the  atmosphere. 


2)  Locally  exponential  atmosphere:  ^  is  constant  over  a 
small  altitude  window. 

3)  P  r  -  constant  atmosphere:  The  dimensionless  quantity 
P  r  remains  constant. 

4)  Isothermal  atmosphere:  The  temperature  is  considered 
constant  through  an  altitude  interval  so  that  dT/dr  =  0  and 
p  becomes  gM/R*T.  For  an  inverse-square  gravitational  field 

g  =  gg  [r^/r]^  (28) 

2 

which  leads  to  the  quantity  P  r  being  constant. 

All  four  profiles  use  the  density  exponential  function  of 
Eq  (22). 

The  biggest  limitation  in  using  this  atmospheric  model  is 
the  assumption  of  spherical  symmetry.  As  expressed  now  the 
air  density  is  only  a  function  of  the  radial  distance  from 
the  center  of  the  Earth.  However,  the  Earth  is  an  oblate 
spheroid  and  hence  produces  a  latitude  effect  to  the 
atmospheric  properties.  This  atmospheric  oblateness  and 
aforementioned  solar  activity  cause  the  main  distortions  from 
spherical  symmetry.  More  complex  models  of  the  atmosphere 
have  been  devised  to  account  for  these  variational  effects. 
See  Reference  4. 


Satel life  Velocity .  The  component  V  in  the  atmospheric 
drag  equation  represents  the  satellite’s  velocity  relative  to 
the  atmosphere.  With  Va  being  the  velocity  vector  of  the 
air  and  v  the  satellite’s  velocity  vector  (both  with  repect 


to  an  Earth  fixed  coordinate  frame) 

V  =  V  -  Va  (29) 

Squaring  Eq  (29)  yields 

=  v^+  Va^-  2  V  Va  cos  y  (30) 

where  y  is  the  angle  between  Va  and  v. 

Let  (J  represent  the  atmospheric  rotation  rate.  This  ang¬ 
ular  velocity  is  assumed  to  be  uniform  about  Earth’s  North- 
South  axis  (27:275).  Thus 

Va  =  r  w  cos  <P  (31) 

where 

r  =  distance  from  Earth’s  center 
<t>  =  geocentric  latitude 

For  small  eccentricities,  a  satellite  travels  nearly  hori¬ 
zontal  through  the  atmosphere  (and  exactly  horizontal  for 
circular  orbits).  A  very  small  error  (<  1%  )  results  if  the 
angle  y  is  taken  as  angle  y*,  the  angle  between  Va  and  the 
horizontal  component  v  of  v  (15:23).  See  Figure  4. 

H 

Examining  spherical  triangle  SNL 

cos  y*  cos  <P  =  cos  i  (32) 

Eq  (31),  with  y  y*,  becomes 

Va  cos  y  =  r  u)  cos  i  (33) 

Substituting  this  result  into  Eq  (30)  produces 

V^  =  v^  f  1  -  (rw/v)  cos  i)^  +  r^ii)^(cos^0  -  cos^i)  (34) 

Further  assumptions  to  Eq  (34)  can  be  made  due  to  the 

variable  and  unknown  rotation  of  the  Earth’s  atmosphere. 

2  2  2  2  2 
First,  the  r  w  term  can  be  neglected  due  to  r  u  <  .005  V 


The  quantity  F  represents  the  effect  of  atmospheric 
rotation  on  the  drag  force.  Assuming  (i)  is  equal  to  Earth’s 
rotation  rate,  F  will  normally  lie  between  0.9  and  1.1  and, 
although  the  effect  of  atmospheric  rotation  on  drag  is 
slight,  it  is  not  inconsequential  (15:25). 

^effici^n^t  oJ|  Drag .  The  drag  coefficient  Cd  is  an 
important  element  in  the  aerodynamic  drag  equation.  While  a 
more  accurate  value  of  the  total  aerodynamic  force  can  be 
determined  using  a  differential  force  equation  it  is  more 
convenient  to  use  the  drag  coefficient  in  the  orbital 
analysis  (19:9).  When  working  with  a  coefficient  of  drag 
several  assumptions  must  be  made  concerning  the  atmospheric 
molecules  (15:14-15): 

1)  The  satellite  is  considered  to  be  stationary  with 
the  air  molecules  flowing  past. 

2)  The  molecules  are  assumed  to  impinge  on  the  satel¬ 
lite,  be  retained  temporarily  on  its  surface,  and  then 
re-em i tted . 

31  The  collisions  between  incident  and  re-emitted 
molecules  are  neglected. 

Several  factors  come  into  play  in  calculating  the  drag 
coefficient  .  The  first  parameter  to  consider  is  the  flow 
regime  through  which  the  satellite  moves.  This  flow  type  is 
determined  by  the  Knudsen  number,  defined  as  the  ratio  of 
mean  free  path  of  atmospheric'  molecules  to  the  charac'ter  i  st  ic 


linear  dimensic^n  of  the  satellite  (11:184).  Two  hundred 


kilometers  above  the  Earth  the  ordinary  continuum  flow  of 
conventional  aerodynamics  no  longer  applies  because  of  the 
extremely  low  air  density  (15:14).  This  region  is  call  free 
molecular  flow  and  has  a  Knudsen  number  of  10  or  greater 
(11:184).  When  a  satellite  is  in  this  flow  regime  its  drag 
coefficient  is  dependent  on  the  molecular  speed  ratio.  This 
is  the  ratio  of  satellite  speed  to  most  probable  molecular 
speed.  For  altitudes  below  700  km  this  speed  ratio  always 
exceeds  5.  This  implies  that  the  random  thermal  motion  of 
the  atmospheric  molecules  can  be  ignored;  i.e.,  the  flow  is 
hyperthermal  (3:931). 

Another  factor  to  the  drag  force  is  the  mechanism  of 
molecular  reflection.  The  energy  exchange  between  the 
atmosphere  molecules  and  the  satellite  is  dependent  on  both 
the  direction  of  the  reflected  molecules  and  their  speed.  It 
is  assumed  the  air  molecules  that  impinge  on  the  satellite’s 
surface  do  not  reflect  specularly  but  instead  attach 
themselves  to  the  outer  layer  of  the  surface  for  a  period  of 
time  before  being  re-emitted.  During  this  period  the 
molecules  ’forget’  their  original  direction  of  motion  and  are 
re-emitted  diffusely.  This  diffuse  reflection  is  strongly 
contingent  upon  the  nature  of  the  satellite’s  surface  and  its 
structure  (3:931). 

The  speed  of  the  re-emitted  molecules  is  determined  by 
their  kinetic  temperature.  During  the  period  of  attachment 
the  molecules  also  ’forget’  their  original  temperatures.  By 
how  much  they  ’forget’  is  uncertain  (15:15)  and  leads  to 


the  accommodation  coefficient,  defined  as 


Ti  -  Tr 
®  “  Ti  -  Te 

where 

Ti  =  original  molecular  temperature 
Tr  =  re-emitted  molecular  temperature 
Ts  =  satellite  surface  temperature 


(38) 


Theoretical  values  of  the  accomodation  coefficient  are  very 
difficult  to  determine.  Many  assumptions  must  be  made 
concerning  the  gas  molecules  and  their  interaction  with  the 
satellite’s  outer  surface.  It  is  suggested  that  low  values 
of  the  coefficient  are  appropriate  (3:931,934).  This  implies 
higher  drag  coefficients.  Reference  15,  on  the  other  hand, 
assumes  from  "conflicting  and  rather  unsatisfactory"  exper¬ 
imental  results  that  the  accommodation  coefficient  is  nearly 
1.0  but  admits  this  assumption  may  be  wrong  (15:15). 

Reference  3  presents  a  graph  of  drag  coefficients  versus 
accommodation  coefficients  for  a  circular  cylinder  with  its 
axis  perpendicular  to  the  direction  of  motion.  Hyperthermal 
free  molecular  flow  and  a  Ts/Ti  ratio  of  0.006  are  both 
assumed  (3:939).  See  Figure  5.  The  upper  line  is  a  olot  of 
Cd  =  2(1+  Ti/6  r)  where  diffuse  re-emission  is  assumed.  The 
term  r  here  is  a  ratio  of  the  speed  of  a  re-emitted  molecule 
\r  to  the  speed  of  an  incident  molecule  Vi  and  is  related 
to  a  by 

r  =  =  [1  +  a  (Te/T.  -1)]’'^  (39) 

The  ratio  Te/Ti  is  very  small  so  r  can  be  approximated  by 


Satellite  Area .  The  projected  area  normal  to  the  free 

molecular  flow  of  Earth’s  upper  atmosphere  also  affects  the 

drag  force  experienced  by  the  satellite.  The  size  of  this 

area  depends  upon  the  orientation  of  the  satellite  as  it 

orbits  the  Earth.  For  a  cylindrical  shaped  body  in  a  gravity 

gradient  stabilized  position  the  projected  area  is  simply 

1  X  d,  where  1  is  the  length  and  d  is  the  diameter.  At  the 

other  extreme,  for  a  cylinder  pointed  into  the  air  flow,  the 

2 

projected  area  is  rr  d  /4. 

However,  a  satellite  might  not  remain  in  these  constant 
area  positions  relative  to  the  flow.  The  system’s  orbital 
dynamics  can  also  alter  the  size  of  the  projected  area. 
Reference  29,  in  a  comparison  of  space  environment  torques, 
concludes  that  below  200  nm  (<"  370  km)  aerodynamic  torques 
dominate  a  satellite’s  dynamics  whereas  above  300  nm  ( -^  555 
km)  the  gravity  gradient  torques  are  the  major  influences. 
Between  these  altitudes  both  space  torques  can  interact  to 
cause  an  unstable  disturbance  in  the  satellite’s  orientation 
(18:2200).  An  uncontrolled  satellite,  with  any  initial 
rotational  motion,  will  soon  begin  spinning  about  its  axis  of 
greatest  moment  of  inertia.  This  action  is  caused  by  these 
relatively  small  external  torques  (15:16). 

With  its  maximum  moment  of  inertia  axis  being  its 
transverse  axis,  the  cylinder  (with  1/d  >  2)  could  assume  an 
orientation  where  its  spin  axis  is  aligned  with  the  air 
flow  (15:16).  This  would  be  its  "aeroplane  propeller"  motion 
and  again  the  projected  area  is  1  x  d.  Conversely,  the 


cylinder  could  have  its  spin  axis  pointing  90  degrees  to  the 

air  flow,  tumbling  end-over-end.  The  mean  cross-sectional 
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area  in  this  case  is  2/7i  (Id  +  1/4  n  d  ).  Of  course  other 
spin  axis  directions  are  also  possible  and  the  projected  area 
would  lie  between  these  two  extreme  values.  Reference  15 
presents  an  expression  for  the  mean  of  these  values  as 

A  =  Id  (0.818+0. 25  d/1)  (42) 

The  two  perturbations  discussed  in  this  chapter  are  the 
major  disturbances  affecting  a  satellite's  orbit  about  the 
Earth.  Each  has  its  own  unique  impact  but  cannot  be 
considered  separately  in  an  orbital  analysis  of  an  External 
Tank.  The  results  of  this  ET  study  will  show  the  coupling 
effect  between  them  that  alters  the  satellite's  orbital  decay 


IV. 


^ o g r mn  I mpl ementat i on 

This  chapter  outlines  the  approach  used  by  the 
Artificial  Satellite  Analysis  Program  to  model  the  Earth’s 
gravitational  field  and  atmosphere.  The  constants  and 
variables  that  describe  the  primary  and  its  gravitational 
field,  the  atmosphere,  and  the  External  Tank  and  its  orbit 
are  also  included. 

Earth,  as  the  primary  in  this  investigation,  is  defined 
by  the  following  parameters  (16:11-14): 


1 ) 

P 

:  3.986006  x  10^ 

2 

km  /sec 

2  ) 

R 

e 

:  6378.14  km 

3  ) 

cj  ; 

e 

:  0.4178074216  x 

_2 

10  deg/sec 

4  ) 

e 

:  0.08199 

e 


ASAP  incorporates  Earth’s  grav'i  tat  i  onal  field  through 
use  of  the  geopotential  expansion  series  Eq  (15).  The  ET  ’  s 
position  in  its  orbit  is  correlated  with  the  associated 
radius,  latitude  and  longitude  of  the  geopotential  and  the 
disturbing  acceleration  is  calculated.  The  values  for  the 
spherical  harmonic  coefficients  C  and  S  of  Eq  (15)  that 
define  Earth’s  mass  distribution  come  from  Table  I. 

For  perturbations  caused  by  a  planet’s  atmosphere  ASAP 
uses  a  simple  exponential  density  model  as  expressed  by  Eq 
(22).  Reference  20  is  the  source  of  Earth’s  atmospheric 
density  and  scale  height  valvaes  at  various  reference 


altitudes.  This  is  a  static  atmosphere  model  (no  time 
dependence)  that  attempts  to  describe  ’average’  atmospheric 
properties  which  represent  a  mean  value  of  the  diurnal, 
seasonal  and  other  solar-caused  variations  detailed  in  the 
previous  chapter.  ASAP  also  assumes  the  atmosphere  is 
rotating  at  the  same  rate  the  Earth  itself  rotates  and 
includes  the  atmos-pheric  oblateness  in  its  calculations. 

The  External  Tank  is  the  main  structural  element  of  the 
Space  Transportation  System,  supporting  both  the  Orbiter  and 
the  two  Solid  Rocket  Boosters  (Gee  Figure  6).  The  tank 
consists  of  three  main  components;  the  LO^  tank,  the 
intertank  and  the  LH  tank.  These  ET  components  are 
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Figure  6.  STS  External  Tank 


31 


manufactured  from  an  aluminum  alloy  and,  after  assembly,  the 


unit  is  coated  with  a  1  to  1-1/2  inch  layer  of  spray-on  foam 
insulation  (8:11-1).  The  following  External  Tank  numbers  are 
used  for  the  orbital  analysis: 


1 ) 

Weight 

:  33503  kg 

2  ) 

Area  (normal  to  flow) 

:  353  m^ 

3  ) 

Cd 

:  2.4 

The  tank  weight  figure  comes  from  Reference  10  and 
represents  an  inert  weight.  It  is  assumed  that  most  of  the 
residual  fuel  at  normal  jettison  altitude  would  be  used  to 
boost  the  Orbiter  and  ET  into  a  higher  orbit. 

According  to  Reference  8,  the  ET,  once  left  in  orbit, 

will  assume  a  gravity  gradient  orientation  where  the 

longitudinal  axis  of  the  tank  points  toward  the  center  of  the 

Earth.  Reference  17  states  that  the  gravity  gradient  torque 

on  an  External  Tank  orbiting  between  400  and  500  km  would 

dominate  all  other  environmental  torques,  including  moments 

caused  by  aerodynamic  forces.  This  situation  will  expose  the 

greatest  amount  of  tank  area  to  the  oncoming  air  molecules. 

Of  course  in  reality,  the  Shuttle  would  jettison  the  ET  with 

its  nose  pointed  into  the  atmospheric  flow  and  there  would  be 

a  certain  amount  of  time  before  the  tank  stabilized  into  the 

gravity  gradient  position.  In  this  case  the  mean  projected 

2 

area  normal  to  the  flow  would  be  less  than  353  m  resulting 
in  less  altitude  loss  over  the  90-day  period.  Howe\er,  for 


the  purposes  of  this  thesis,  the  worst  case  for  the  projected 


area  will  be  considered  where  the  External  Tank  is  dropped 
off  directly  into  a  gravity  gradient  position.  This  ET  area 
is  calculated  from  dimensions  given  in  Reference  26. 

For  cylindrical  shaped  satellites,  drag  coefficients 
derived  in  References  3,9,11,15  and  22  range  in  value  from 
2.0  to  3.0.  Modeling  the  External  Tank  as  a  cylinder,  a  Cd  of 
2.4  is  used  as  a  representative  figure. 

The  classical  Keplerian  elements  (  a  ,  e  ,  i  ,  0, o),  M  )  are  used 
to  define  the  initial  starting  conditions  of  the  orbiting 
External  Tank.  During  preliminary  orbit  investigations 
various  values  of  and  M  were  evaluated  and  found  to  have 

no  discernable  effect  on  the  90-day  altitude  loss  rate.  For 
this  reason,  during  all  subsequent  test  runs,  these  three 
elements  are  set  to  zero  for  the  initial  orbit  condition. 


V.  Results 


The  results  of  this  thesis  are  discussed  in  two 
sections.  The  first  section  details  the  search  for  the 
lowest  possible  altitude  for  the  defined  External  Tank 
parking  orbit:  Maximum  25  km  altitude  loss  over  90  days. 

The  orbital  elements  of  this  minimum  altitude  orbit  are  then 
explored.  The  second  section  investigates  various  aspects  of 
the  orbit,  gravitational  field,  atmosphere  and  ET  character¬ 
istics  to  determine  their  effect  on  the  magnitude  of  altitude 
loss.  An  introductory  section  is  needed  however  to  define 
’altitude  loss’  as  used  in  this  thesis. 


Altitude  Loss  Definition 


Figure  7  is  a  plot  of  the  daily  fluctuations  in 
geocentric  altitude  of  an  orbiting  External  Tank  perturbed  by 
both  a  gravity  field  and  atmospheric  drag.  These  fluctuations 
do  not  lead  to  an  easy  statement  of  the  amount  of  altitude 
lost  by  the  ET  over  the  90-day  period.  However,  a  regression 
function  is  used  to  perform  a  linear  least  squares  fit  to  the 
data.  See  Figure  8.  The  program  calculates  the  slope  of 
this  line  which  represents  the  mean  altitude  loss  rate.  It 
is  therefore  very  convenient  to  define  the  altitude  loss  in 
kilometers  as: 

90-Day  Altitude  Loss  =  Line  Slope  ( km/day )  x  90  (days) 


Figure  7  . 


Geocentric  Altitude  Fluctuations 


Minimum  Park i ng  Orbit  A1 t i tude 

To  perform  this  phase  of  the  External  Tank  study  initial 

O 

conditions  of  28  inclination  and  zero  eccentricity  were  used 
A  3x3  gravitational  field  and  a  simple  atmosphere  model  were 
integrated  into  this  analysis.  The  search  for  the  minimum 
altitude  began  at  500  km  and  progressed  downward  at  25  km 
increments.  Computer  test  runs  were  conducted  until  the 
parking  orbit  constraint  had  been  violated.  Results  of  this 
investigation  are  reported  in  Table  II.  Due  to  the  exponen- 


REGRESSION  OF  GEOCENTRIC  ALTITUDE  DATA 


Figure  8.  Linear  Regression  to  Geocentric  Altitude  Data 

tial  nature  of  the  atmospheric  density,  there  is  a  substan¬ 
tial  increase  in  the  altitude  loss  rate  as  the  initial 
altitude  of  the  ET  orbit  decreases.  The  altitude  of  425  km 
meets  the  established  parking  orbit  criteria. 

T^ime  H  i  story  of  the  Orbital  Elements 

Using  425  km  as  the  initial  altitude  for  the  External 
Tank  orbit,  the  orbital  elements  were  investigated  for  a 


period  of  90  days.  Due  to  the  fluctuations  in  the  elements 


TABLE  II 

90-Day  Altitude  Loss  for  Initial  ET  Altitudes 


Initial  Orbit  Altitude  (km) 

500 

475 

450 

425 

400 


90-Day  Altitude  Loss  (km) 

5 . 5 
8.4 

13.2 

21.8 

38 . 2 


a  curve  fitting  program  was  adapted  from  Reference  7.  This 
technique  suppresses  the  short  term  variations  in  the  orbital 
elements  and  allows  the  long  period  behavior  to  be  seen. 

The  program  calculates  the  coefficients  to  1st,  2nd  and  3rd 
order  polynomials  and  determines  which  order  equation 
produces  the  best  fit  to  the  data.  These  polynomials 
developed  for  the  orbital  elements  represent  a  convenient 
method  to  calculate  the  External  Tank’s  orbit  at  some  future 
time  t,  within  the  90-day  period  studied.  Figures  9  through 
13  are  plots  of  the  daily  variation  in  the  orbital  parameters 
and  the  associated  long  period  trend  line. 

Semi-Major  Axis  (Figure  9).  The  semi-major  axis  of 
this  External  Tank  parking  orbit  contracts  by  21.4  km 
( 0.3%)  during  the  90-day  period.  The  long  term  trend  line 
is  a  2nd  order  polynomial: 

a  =  6801.2279327  -  0.2122813  t  -  0.0003232  t^  (43) 
Polynomials  generated  for  other  initial  altitudes  show  that  a 


SEMI-MAJOR  AXIS  TIME  HISTORY 

(3x3  Gravity  Flold/Atmospharlc  Drag) 


Figure  9.  Semi-Major  Axis  Time  History 

1st  order  equation  is  adequate  for  500  and  475  km  but  start¬ 
ing  at  450  km  a  2nd  order  polynomial  is  required  for  the 
best  fit  to  the  semi-major  axis  data. 

Eccentricity  (Figure  10),  Although  the  External  Tank 
was  placed  in  a  circular  orbit  (e  =  0)  for  this  investi¬ 
gation,  the  value  of  this  orbit’s  eccentricity  did  not  remain 
at  zero.  This  variation  in  eccentricity  is  produced  by  the 
gravitational  perturbations.  The  best  fit  polynomial  is 


Longitude  of  Ascending  Node  (Figure  12).  The  best  fit 
polynomial  for  this  orbital  element 

ft  =  360.1578150  -  7.0659969  t  (46 

exactly  matches  the  data  generated  by  ASAP.  This  graph  has 
the  node  initialized  at  360  and  indicates  the  External 
Tank’s  orbital  plane  rotates  about  Earth’s  polar  axis  approx 
imately  once  every  50  days  or  about  7°per  day  This  effect  on 
ft  is  caused  by  gravitational  perturbations  (15:7-8). 


ASCENDING  NODE  TIME  HISTORY 

(3x3  Gravtty  F]«ld/Atmoaph«r1c  Drag) 


TIME  (day*) 
0  ASAP  Data 


Figure  12.  Longitude  of  Ascending  Node  Time  Hist 


Q.®  a  1 1  i  tude  (Figure  13).  The  daily  variation  in 

the  geocentric  altitude  of  the  External  Tank  is  the  result  of 
the  varying  force  of  attraction  on  the  tank  caused  by  Earth’s 
gravity  field.  The  best  fit  polynomial  is  1st  order 

Geo.  Alt.  =  417.1935674  -  0.2459279  t  (47) 

and  shows  an  altitude  loss  rate  of  about  0.25  km  per  day. 

Orbital  element  polynomials  for  other  altitudes  are 
contained  in  Appendix  B.  The  next  section  of  the  results 
will  discuss  the  different  influences  on  the  External  Tank 
altitude  loss  rate. 


GEOCENTRIC  ALTITUDE  TIME  HISTORY 

(3x3  Gravity  n.ld/AtTne.ph.r1c  Drag) 
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Figure  13.  Geocentric  Altitiide  Time  Historv 


Initial  Al t i tude  Effects 

Ninety-day  altitude  loss  data  from  the  parking  orbit 
investigation  in  Section  I  was  used  to  produce  Figure  14.  At 
500  km  the  upper  atmospheric  density  is  extremely  low  and  the 
External  Tank  loses  very  little  altitude  during  the  90  days. 
However,  at  lower  altitudes  the  atmospheric  drag  becomes  more 
of  a  factor  and  shrinks  the  ET ' s  orbit  considerably. 


ORBIT  ALTITUDE  EFFECTS  ON  ALTITUDE  LOSS 
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Figure  14.  Initial  Altitude  Effects  on  orbit  I) 
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Orbital  Elejneiit  Ef  f  ects  :  Incl  ination  and  Eccent  r  ic  i  ty 

Both  these  investigations  were  conducted  using  a  3x3 
gravitational  field  and  a  simple  atmosphere  model.  The 
initial  altitude  of  the  External  Tank  orbit  was  425  km. 

For  the  inclination  study  test  runs  were  made  at  18,  23, 
28,  33,  and  38  degrees  with  zero  initial  eccentricity  in  each 
case.  Figure  15  indicates  the  slight  increase  in  the  90-day 
altitude  loss  as  the  orbit  inclination  is  decreased.  This 
effect  is  due  to  the  oblateness  of  the  modelled  atmosphere 
where,  for  a  given  altitude,  the  air  is  slightly  thicker  at 
the  equator  than  at  the  poles. 


INCLINATION  EFFECTS  ON  ALTITUDE  LOSS 
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Figure  15. 


Inclination  Effects  on  Orbit  Decav 


For  the  next  study  test  cases  were  run  with  eccentricity 
values  of  0.0,  0.005,  0.0075  and  0.01  with  an  initial 

O 

inclination  of  28  in  each  case.  As  seen  from  Figure  16,  by 
increasing  the  initial  eccentricity  of  the  ET’s  orbit,  the 
90-day  altitude  loss  increases.  The  explanation  of  this 
effect  is  found  by  looking  at  the  variation  in  the  geocentric 
altitude  which,  in  turn,  determines  the  atmospheric  density. 
For  near  zero  eccentricity  orbits,  the  value  of  geocentric 
altitude  will  vary  as  was  shown  in  the  first  section  of  the 
results.  However,  as  the  eccentricity  is  increased,  it  is 


ECCENTRICITY  EFFECTS  ON  ALTITUDE  LOSS 
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Figure  16.  Eccentricity  Effects  on  Orbit  Decay 
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found  that  the  fluctuations  in  geocentric  altitude  become 


larger.  This  results  in  the  External  Tank  being  exposed  to 
higher  values  of  air  density  in  addition  to  having  a  higher 
orbital  velocity  at  the  lower  extremes.  These  two  factors 
result  in  a  greater  drag  force  as  expressed  by  Eq  (20). 


Gravitational  Field  Effects 


In  the  first  part  of  this  study  the  External  Tank  was 
put  in  an  orbit  with  initial  cond itions  of  28  inclination 
and  zero  eccentricity  with  computer  runs  made  at  initial 
altitudes  of  425  and  500  km.  For  each  altitude  90-day  orbits 
were  calculated  using  a  simple  atmosphere  model  with  and 
without  a  3x3  gra\ity  field.  Table  II  details  the  coupling 
effect  between  these  two  orbit  perturbations.  Due  to  the 
varying  gravitational  attraction  of  a  3x3  field,  the  ET  finds 
itself  at  times  pulled  into  a  lower  than  normal  orbit  thereby 
experiencing  greater  atmospheric  drag  forces.  The  overall 
effect  of  a  complex  geopotential,  combined  with  an  atmos¬ 
phere,  is  a  higher  rate  of  altitude  loss. 


TABLE  III 


Air  Drag/Gravity  Field  Coupling  Effects 


A 1 1 i tude  ( km ) 


90-Dav  Altitude  Loss  (km) 


A i r  Drag 


Air  Drag  +  3x3  Grav  Field 


25 


The  next  investigation  into  the  gravitational  field 
effects  on  altitude  loss  concerns  the  degree  of  complexity 
needed  for  acceptable  results.  Using  an  initial  External 

o 

Tank  orbit  of  425  km  altitude,  28  inclination  and  zero 
eccentricity,  test  cases  were  run  with  a  simple  atmosphere 
and  models  of  the  following  geopotential  fields:  2x0,  2x2, 
3x1,  3x3.  In  all  four  cases  the  90-day  altitude  loss  was 
21.8  km.  The  2x0  field  models  the  oblate  shape  of  the  Earth 
and  is  the  dominant  gravitational  perturbation. 


At mospher ic  Dens i ty  Effects 


The  simple  exponential  atmosphere  model  is  one  of  the 
major  assumptions  made  in  this  orbital  analysis  of  an 
External  Tank.  This  static  model  leaves  out  the  dynamic 
characteristics  of  Earth’s  atmosphere  as  discussed  in  Chapter 
III.  All  the  fluctuations  in  the  atmospheric  properties  are 
averaged  to  produce  a  mean  density.  However,  over  any 
designated  time  period,  the  value  of  this  mean  density  at  any 
given  altitude  could  be  much  higher  than  that  sited  by 
Reference  20,  Subsequently,  a  series  of  computer  runs  were 
done  to  study  the  effects  of  increased  atmospheric  density  on 
an  External  Tank’s  orbit.  Density  numbers,  with  associated 
scale  height  values,  were  taken  from  Reference  12  covering  a 
low-density  to  a  medium-density  atmosphere.  An  ET  orbit  of 

o 

425  km  altitude  ,  28  inclination  and  zero  eccentricity 


m 


formed  the  basis  for  studying  the  90-day  altitude  loss.  As 


was  the  case  when  the  various  starting  altitudes  were 
examined,  a  higher  mean  value  of  air  density  will  cause  the 
ET  ’  s  altitude  to  decre^ise  markedly.  See  Figure  17. 


DENSITY  EFFECTS  ON  ALTITUDE  LOSS 

(425km  OrbIt/SxS  Gray  Fisid/Atmei  Drag) 

70  -i 


0  4- - 1 - 1 - 1 - 1 - r 

2.09E-04  9.41  E-04  1.67E-03  2.40E-03  3.14E-03  3.87E-03  4.60E-03 

ATMOSPHERIC  DENSITY  (kg/cubIc  km) 


Figure  17.  Atmospheric  Density  Effects  on  Orbit  Decay 


External  Tank  Characteristics  Effej:ts 


The  values  used  in  this  thesis  for  the  External  Tank’s 
projected  area  to  the  air  flow  and  it's  drag  coefficient  are 
valid  assumptions  but  by  no  means  certainties.  As  was 


dis(,'ussed  in  Chapter  III  a  combination  of  gra\’ i  tai  onal  and 


aerodynamic  forces  can  produce  different  orientaions  of  an 
orbiting  spacecraft  which,  in  turn,  will  affect  the  magnitude 


of  the  projected  area  and  drag  coefficient.  With  this  in 
mind  a  study  was  done  using  the  ballistic  coefficient  as  a 
parameter  to  vary.  From  Eq  (38) 

Cd  A 


B  = 


2m 


{  48) 


For  an  ET  kept  with  its  front  end  into  the  air  flow,  the 


projected  area  would  be  55  m  .  Conversely,  the  gravity 

2 


gradient  profile  has  an  approximate  area  of  353  m 
References  23  and  30  give  higher  Cd  values  for  various 
orientations  of  a  cylindrical  satellite  so  a  Cd  range  of  2  to 
5  is  used,  along  with  an  ET  mass  of  33503  kg,  to  calculate  a 
minimum  and  maximum  ballistic  coefficient.  The  minimum  would 
be 


Bm  i  n  = 


(2)  (55) 


(2)  (33503) 


=  1.64  X  10'^  m^/kg 


and  the  maximum 


„  _  (  5)  (353)  _  „  ,„-2  _2,,  ^ 

=  (2)  ■(  33503)  =  ^ 

Computer  runs  were  done  using  ballistic  coefficients  of 
l,64xl0;^  9.86xl0;^  1.81xl0*^and  2.63x10'^  m^/kg  convolved 
with  initial  altitudes  of  400,  425,  450  and  475  km,  and  with 

O 

nominal  values  of  28  inclination  and  zero  eccentricity. 

This  convolution  of  ballistic  coefficients  and  starting 
altitudes  produced  16  values  of  90-day  altitude  loss.  To 
develop  an  idea  of  the  External  Tank  characteristics  and 


their  effect  on  the  loss  rate  a  3-dimensional  plot  was 


formed  the  basis  for  studying  the  90-day  altitude  loss.  As 
was  the  case  when  the  various  starting  altitudes  were 
examined,  a  higher  mean  value  of  air  density  will  cause  the 
ET ’ s  altitude  to  decrease  markedly.  See  Figure  17. 
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Figure  17.  Atmospheric  Density  Effects  on  Orbit  Decay 


External  Tank  Characteristics  Effects 


The  values  used  in  this  thesis  for  the  External  Tank's 
projected  area  to  the  air  flow  and  it’s  drag  coefficient  are 
valid  assumptions  but  by  no  means  certainties.  As  was 
discussed  in  Chapter  III  a  combination  of  gravitaional  and 
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accomplished  (See  Figure  18.  Note:  Due  to  the  graphics 
program  operation  the  ballistic  coefficients  labelled  on  this 
plot  should  be  multiplied  by  2x10  ^). 


Figure  18.  3-D  Graph  of  90-Day  Altitude  Loss 

Figure  19  is  a  contour  plot  of  the  same  data  giving  a 
clearer  picture  of  the  90-day  altitude  loss  for  various 
combinations  of  ballistic  coefficient  and  initial  External 
Tank  orbit  altitude. 
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Figure  19.  Contour  Graph  of  90-Day  Altitude  Loss 

Th i rd  Body  Effects 

Since  the  Artifical  Satellite  Analysis  Program  has  the 
capability  to  include  third  body  perturbations,  an  investi¬ 
gation  was  conducted  to  determine  if  the  Sun  or  Moon  would 
have  any  effect  on  an  orbiting  External  Tank.  Teat  cases  were 


performed  with  the  standard  orbit  of  28°  inclination,  zero 
eccentricity  and  an  initial  altitudes  of  425  km.  Both 


Earth’s  gravitational  and  atmospheric  perturbations  were 
included.  The  third  body  situations  explored  were:  Sun  only 
Moon  only  and  Sun/Moon  together.  Table  IV  lists  the  curve 
fitting  polynomials  calculated  for  the  Sun/Moon  computer  run 
No  noticeable  differences  were  found  between  the  orbit 
histories  produced  here  and  the  like  profile  accomplished  in 
the  first  section  of  the  results. 

TABLE  IV 

Orbit  Element  Polynomials  for  Sun/Moon  Effects 

a  =  6801.2284434  -  0.2122904  t  -  0.0003230  t^ 

e  =  0.0012659  +  0.0000010  t 

i  =  27.9840822  -  0.0000204  t 

0  =  360.1587464  -  7.0663866  t 

Geo. Alt.  =  417.1945712  -  0.2459292  t 


VI  . 


Conclusions  and  Recommendations 


Conclusions 

This  orbital  analysis  of  an  External  Tank  in  low  Earth 
orbit  represents  a  'first  look’  at  how  the  tank  would  react 
to  the  perturbative  effects  of  Earth’s  gravitational  field 
and  atmosphere.  A  scenario  involving  an  ET  parking  orbit  was 
developed  as  a  way  to  limit  the  scope  of  the  problem.  The 
gravity  gradient  attitude,  the  coefficient  of  drag  and  the 
simple  atmosphere  model  were  the  major  assumptions  to  this 
investigation . 

A  major  outcome  of  this  study  was  the  425  km  altitude, 
determined  to  be  the  lowest  initial  altitude  for  keeping  the 
External  Tank  in  the  specified  parking  orbit.  Under  the 
presumptions  used  in  the  scenario,  the  ET  lost  approximately 
22  km  during  the  90-day  period  examined.  The  425  km  altitude 
is  within  reach  of  the  Shuttle/ET  combination  (5:2),  thus  no 
ET  reboost  would  be  necessary  under  these  stated  conditions. 

Due  to  the  daily  variations  in  the  calculated  orbital 
elements  a  curve  fitting  program  was  used  to  extract  their 
long  term  behavior.  In  this  process  the  Artificial  Satellite 
Analysis  Program  generated  the  time  histories  of  the  orbital 
elements  and  an  ’averaging*  was  then  performed.  The 
polynomials  produced  by  using  ASAP  output  data  provide  a 
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quick  method  of  determining  future  values  of  the  External 
Tank  orbital  elements. 

From  this  initial  study  of  an  External  Tank  in  space, 
several  conclusions  can  be  drawn  from  the  long  period  trends 
in  the  ET ’ s  orbital  parameters.  As  with  all  satellites  in 
LEO,  the  atmosphere  causes  the  semi-major  axis  to  shrink  at 
an  ever  increasing  rate.  When  deposited  in  orbit  with  zero 
eccentricity,  the  ET  will  remain  in  a  near  circular  orbit. 

The  inclination  of  its  orbit  will  also  stay  close  to  initial 
conditions.  The  orientation  of  the  ET  orbital  plane,  as 
defined  by  the  longitude  of  ascending  node,  cycles  approx¬ 
imately  every  50  days.  Any  rendezvous  missions  with  an 
on-orbit  External  Tank  will  have  to  take  this  into  account. 

The  examination  of  the  various  factors  that  influence 
the  External  Tank  orbit  decay  rate  showed  that  some  factors 
are  more  detrimental  than  others.  The  ET  inclination  and 
eccentricity  parameters,  within  the  nominal  Shuttle  operating 
limits,  have  minor  effects  on  orbit  contraction.  On  the 
other  hand,  the  3-dimensional  graph  and  associated  contour 
plot  (Figures  18  and  19)  demonstrate  the  combined  effects  of 
atmospheric  density  and  ET  characteristics  of  area,  mass,  and 
drag  coefficient.  At  lower  altitudes  where  the  air  density 
is  greater,  small  increases  in  the  ballistic  coefficient 
produce  large  changes  in  the  altitude  lost  by  the  External 


Tank  during  the  90-day  period. 


i 
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The  Earth’s  gravity  field  cannot  be  neglected  in  an 
analysis  of  an  orbiting  External  Tank.  The  gravitational  and 


atmospheric  drag  perturbations  interact  to  produce  a  greater 


altitude  loss  than  their  individual  effects.  However,  for  a 
long  term  orbital  study,  the  short  period  gravitation 
anomalies  can  be  disregarded  since  a  2x0  field  that  models 
Earth’s  oblateness  can  provide  sufficient  accuracy. 


Recommendations 


As  with  any  initial  study  into  a  problem,  follow-on  work 
is  appropriate.  Several  diverse  areas  could  be  examined  for 
greater  insight  into  an  orbiting  External  Tank.  As  modelled, 
the  atmosphere  used  in  this  thesis  was  one  of  the  simplest 
possible.  It  would  be  beneficial  to  employ  a  dynamic  model 
where  the  time  dependence  of  the  atmospheric  properties  is 
considered.  To  be  on  the  conservative  side  this  study 
started  the  ET  in  a  gravity  gradient  position  which  presented 
the  largest  possible  area  to  the  air  flow.  At  some  point  in 
the  ET’s  orbit  decay,  moments  due  to  aerodynamic  forces  will 
become  the  predominant  torque  on  the  tank  and  it  will  then 
transition  from  a  stabilized  orientation  into  an  uncontol- 
lable  state.  An  investigation  into  the  ET’s  orbit  dynamics 
would  provide  a  more  accurate  value  for  its  projected  area. 
Also,  an  analysis  could  be  accomplished  to  determine  the 
amount  of  lift  generated  by  the  ET  and  the  effect  on  orbital 
contract i on . 
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In  addition  to  this  projected  area  another  uncertainty 
is  the  External  Tank’s  coefficient  of  drag.  The  value  of  Cd 
is  dependent  on  the  accommodation  coefficient  where  the 
interaction  of  the  air  molecules  with  a  satellite’s  surface 
plays  an  important  part.  However,  a  foam  covered  ET  presents 
quite  a  different  exterior  to  the  atmospheric  flow  than  the 
normal  metallic  surfaces  of  other  spacecraft.  An  examination 
into  this  area  could  result  in  different  values  of  the  drag 
coefficient  being  appropriate  for  use  in  an  ET  orbital 
analysis.  As  emphasized  earlier,  even  small  changes  in  these 
External  Tank  characteristics  can  have  considerable  influence 
in  the  orbital  decay  rate. 

The  possibility  exists  that  multiple  External  Tanks 
could  be  stored  in  an  orbiting  ’tank  farm’  (P'igure  20). 


Looking  at  the  drag  deceleration  equation 
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m 


Cd  A  ,,2 
P  V 


2m 


(  21  ) 


the  magnitude  of  the  deceleration  could  be  decreased  by  the 
proper  orientation  of  two  or  more  tanks.  It  may  be  possible 
to  fix  the  attitude  of  the  multiple  tanks  so  the  increase  in 
the  area  exposed  to  the  air  flow  is  more  than  offset  by  the 
increased  tank  mass.  However,  an  analysis  of  multi-tank  drag 
coefficients  would  be  required  first. 

Another  possible  avenue  of  study  would  be  to  dev^elop  the 
method  of  averages  analytical  technique  discussed  in  Chapter 
III  to  include  both  gravitational  and  atmospheric  drag 
perturbations.  A  comparasion  could  then  be  done  between  this 
technique  and  the  approach  followed  by  this  thesis. 


The  proposed  uses  of  an  orbiting  External  Tank  bolster 
the  idea  that  it  is  too  valuable  a  resource  to  be  expended  on 
each  Space  Shuttle  flight.  The  United  States  has  previously 
adapted  other  space  hardware  for  new  applications.  The  prime 
example  of  this  flexibility  is  Skylab,  where  a  third  stage 
from  a  Saturn  V  was  converted  into  a  space  station.  But, 
like  Skylab,  an  External  Tank  will  not  stay  in  orbit 
indefinitely.  Plans  and  decisions  must  be  made  prior  to  its 
use  in  space.  This  thesis  has  attempted  to  provide  some 
preliminary  answers  to  this  problem. 


Appendix  A 


Conversion  of  the  Spherical  Harmonic  to  Kepler ian  Elements 


This  section  will  be  concerned  with  transforming  the 
expression  for  the  geopotential  discussed  in  Chapter  III  into 
a  form  that  can  be  used  in  the  method  of  averages  and  is 
developed  mainly  from  References  2  and  14.  The  potential  is 
a  function  of  radius  r,  latitude  4>  and  longitude  X: 
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For  use  in  the  method  of  averages  this  potential  needs  to 
converted  to  the  Keplerian  elements  a,e,i,Q,u  and  m. 

To  accomplish  this  element  transformation  several 
trigonometric  identities  will  be  required: 
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One  term  of  the  potential  series  can  be  represnted  as 
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S  where  a  is  the  equatorial  radius  of  body  M. 
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Working  with  Figure  21  a  substitution  is  now  made  for  mX, 

mX  =  [m(a-ft)  +  m((>-6)]  (8) 

where  a  is  the  right  ascension  and  6  is  Greenwich  sideral 
time;  i.e.  X  =  a-6.  This  produces 
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Figure  21.  Satellite  Orbital  Elements 


Using  spherical  trigonometric  relationships  from  Figure 
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where  f  is  the  true  anomaly. 

From  Eqs  ( 10 )  and  (11)  come 
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Now;  trigonometric  identities  (4)  are  applied  to  Eqs  (9) 


and  identities  (12)  are  substituted  into  the  result.  These 
steps  yield 
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The  terms  on  the  right  side  of  Eq  (13)  are  substituted 
for  sin<l'  into  the  associated  Legendre  function  Eq  (2).  This 
new  F  with  Eqs  (14)  are  then  placed  into  the  potential  Eq 
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where  k  is  the  integer  pa rt  of  (l-m)/2.  ldentit\  (o)  along 
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Applying  identities  (6)  to  products  of  trigonometric 
functions  in  Eq  (16)  and  neglecting  any  term  with  an  odd 
power  of  j  (since  V  is  real)  results  in 
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where  k  is  the  integer  part  of  (l-m)/2.  Also,  t  is  summed 
from  zero  to  the  lessor  of  p  or  k  and  c'  is  summed  over  all 
values  making  the  binomial  coefficients  nonzero.  'I’hese 
summations  for  t  and  c  are  due  to  the  definition  of  the 
Binomial  Coefficient.  The  lower  term  in  the  coefficient  must 
be  equal  to  or  greater  than  .zero  and  the  upyier  t.erm  must  be 
ei]ual  to  or  greater  th.an  the  lower  term. 

At  this  point  the  potential  equation  still  depends  on 
the  radius  r  and  the  true  anomaly  f.  These  coordinates  ('an 
be  replaced  with  the  semi-ma.jor  axis  a,  ecf'eri  t  r  i  c  i  t  y  f  and 
mean  anomaly  M.  Extracting  the  general  trigonometric  term 
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where  x"'""  are  called  Hansen’s  coefficients.  Substituting 
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Eqs  (22)  (with  t  =  l-2p  +  q,  m  ::  l-2p  and  n  =  -1-1)  into 
expression  (21)  yields 
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From  this,  tho  er^npral  term  (20)  1  potjmes 
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Determining  Hansen’s  coefficients,  x"’"',  is  quite  an  involved 
process.  Their  development  can  be  found  in  Reference  1. 

The  coefficients,  X  *  ,  are  dependent  on  eccentricity 
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alone  so  an  Eccentricity  function  G  (e)  can  be  defined: 
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At  this  point,  with  the  preceding  der'init  ion  lAjr 
Hanser\’s  coe  1' f  i  c  i  ent  s  ,  the  general  term  (  z!  4)  ';an  be  inserted 
back  into  the  potential  equation  (18)  producing 


p  a '  1  t 

V  r  - — ^  ^  F,  (i)  £  G.  (e)  S,  ((J,M,O,0) 

+  1  ^  Imp  I  p  q  1  in  p  ri 


1  m 


1 


p  =  0 


(  Z9  ) 


where 


I  m  p  q 


I  m 


I  m 


1  -  m  even 


cos  [(l-2p)w  +  (l-2p  +  q)M  +  m(0-(4)J 

I  -  m  odd 


-1  I  -m  ev  e  n 


I  m 


I  m 


sin  [(l-2p)g)  +  (l-2p+q)M  +  m(Q-O)]  (JO) 


1  -  m  odd 


A  table  of  F  (i)  and  G  (e)  functions  is  given  in 
I mp  1 pq 

Reference  14.  This  completes  the  transformation  of  the 
potential  from  spherical  to  Keplerian  coordinates.  At  this 
stage,  the  geopotential  V  c'ould  be  a\'eraged  over  one  orbital 
period  (12:3-4): 

1  r‘'” 

<  Vav^-  >  r  -  V  ( M )  dM  (3  1) 

Zn 

This  would  eliminate  any  depi'ndenc'f'  of  Hq  (  1  )  on  the  mean 
anomaly  and  remove  all  short-term  variat ions  from  this 
disturbing  function  and  subsequently,  from  the  equations  of 


mot  i  ( )  n . 


Appejidix  B 


Orbital  Element  Polvnomialb 

The  following  polynomials  Vvere  generated  from  External 
Tank  orbit  data  c'a  it'u  1  at.ed  by  the  Ai-tificial  Satellite 
Analysis  Pr'cigram  at,  i  r>  i  t  i  a  1  altitudes  of;  400,  4Z5,  450,  475 

and  500  kilometers. 


4  00  km 

a  =  57  75.9:568930  -  0.3960781  t  -  0.0014  187  t" 

e  =  0.0013940  -  0.0000023  t 

1  =  27.9842948  -  0.0000305  t 

0  :  360,2713382  -  7.1702218  t 

Geo. Alt  =  391.4884943  -  0.3029662  t  -  0.0013949  t" 


H  =  6801.2279327  -  0.2122813  t  -  0.0003232  t" 

e  =  0.0012661  +  0.0000010  t 

i  =  27.9844744  -  0.0000193  t 

0  -  360.1578150  -  7.0659969  t 

Geo.  Alt  =  4  17.1935674  -  0.2159279  t. 


s==:^^=======d:^^=^:£===s==z====: 

a  =  6825.9871108  -  0.1213986  t  -0.0003565  t“ 

e  :  0.0014040  -  0.0000010  t 

i  =  27.9858543  -  0.0000615  t 

0  -  360.1058726  -  6.9688603  t 

Geo. Alt.  =  441.8590712  -  0.1519032  t 


6  ( 


a 


6851 .216971 1 


0,0949190  t 


e  =  0.00139  32  -  0.0000004  t, 

i  =  27,9846431  -  0.0000141  t 

0  r  360.0663384  -  6,8757352  t 
Gf-o.Alt.  r  466.7700440  -  0.0982593  t 


5  00  km 

a  =  6876.1566632  -  0.0619361  t 

e  =  0.0013675  +  0.0000002  t 

i  =  27.9848512  -  0.0000140  t 

0  =  360,0416657  -  6.7857346  t 

Geo. Alt  =  491.7578786  -  0.0658334  t 
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Abstrac t 


By  flying  a  different  launch  profile  ,  it  is  possible 
for  the  Space  Transportation  System’s  Orbiter  to  bring  the 
External  Tank  directly  into  space.  Many  studies  by  NASA  and 
private  industry  have  detailed  the  potential  on-orbit  uses  of 
an  External  Tank.  However,  at  Space  Shuttle  operating 
altitudes,  an  orbiting  tank  will  experience  multiple 
environmental  forces  resulting  in  its  decay  into  the  lower 
atmosphere  and  eventual  re-entry. 

This  thesis  conducts  a  preliminary  study  of  a  single 
External  Tank  in  low  Earth  orbit.  Criteria  for  a  parking 
orbit  are  defined  and,  using  an  orbit  prediction  computer 
program  with  atmospheric  drag  and  gravitational  perturbations 
included,  a  search  is  made  for  the  lowest  initial  altitude 
that  will  allow  the  External  Tank  to  remain  in  this  orbit 
window.  The  starting  altitude  that  meets  the  orbit 
requirements  is  found  to  be  within  reach  of  the  Shuttle’s 
capabilities.  The  orbital  elements  of  this  parking  orbit 
are  then  analysed  and  a  method  for  quick  calculation  of  these 
fiarameter’s  is  demised.  An  evaluation  of  the  factors  that 
affect  the  orbital  contraction  of  an  External  Tank  is  also 
performed.  The  atmospheric  density  and  the  tank 
<:ha  rac  t  e  r  1  3  1 1  cs  can  both  contribute  to  high  orbital  decay 
rat  t  *  s  . 


